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ABSTRACT 

We have looked for bulk motions of galaxy clusters in the WMAP 7 year data. We isolate the 
kinetic Sunyaev-Zeldovich (SZ) signal by filtering the WMAP Q, V and W band maps with multi- 
frequency matched filters, that utilize the spatial properties of the kinetic SZ signal to optimize 
detection. We try two filters: a filter that has no spectral dependence, and a filter that utilizes 
the spectral properties of the kinetic and thermal SZ signals to remove the thermal SZ bias. We 
measure the monopole and dipole spherical harmonic coefficients of the kinetic SZ signal, as well 
as the £ = 2 — 5 modes, at the locations of 736 ROSAT observed galaxy clusters. We find no 
significant power in the kinetic SZ signal at these multipoles with either filter, consistent with 
the ACDM prediction. Our limits are a factor of ~ 3 more sensitive than the claimed bulk flow 
detection of Kashlinsky et al. (2009). Using simulations we estimate that in maps filtered by our 
matched filter with no spectral dependence there is a thermal SZ dipole that would be mistakenly 
measured as a bulk motion of ~ 2000 — 4000 km/s. For the WMAP data the signal to noise ratio 
obtained with the unbiased filter is almost an order of magnitude lower. 

Subject headings: Cosmology: cosmic background radiation, observations, diffuse radiation 



1. Introduction 

Verifying the depth of convergence of the 627 
km/s motion of the Local Group (LG) with respect 
to the frame defined by the Cosmic Microwave 
Background (CMB) has been a long-standing goal 
of observational cosmology. Cosmographers seek 
to identify the region in which the local universe 
is at rest with respect to the CMB frame, and the 
LG motion has converged to the CMB dipole di- 
rection. This region would contain all of the mass 
sources responsible for the motion of the LG and 
would define a sample of the universe representa- 
tive of the whole. The coherent motion of matter 
caused by gravitational potentials is termed a bulk 
flow (e.g., Strauss & Willick 1995); if the flow cov- 
ers the entire universe it is sometimes referred to 
as a dark flow (e.g., Kashlinsky et al. 2008). Such 
a flow is equivalent to an intrinsic dipole of the 
universe. Models of inflation can be constructed 
to explain such a measurement (e.g, Turner 1991; 
Kashlinsky et al. 1994). 
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To date there is no consensus on the depth of 
convergence. The CMB dipole magnitude and di- 
rection have been precisely measured by WMAP: 
(3.355 ± 0.008) mK in the direction of galactic 
longitude I = 263.99° ± 0.14° and latitude b = 
48.26° ± 0.03 (Jarosik et al. 2010). The LG mo- 
tion towards the Great Attractor was measured 
by Dressier et al. (1987) in the direction of / = 
312°, 6 = 6° and found to be coherent out to 30ft -1 
Mpc and near zero by 40/i _1 Mpc. Using a full- 
sky peculiar velocity survey, with a depth ranging 
from 80/i _1 Mpc to llO/i" 1 Mpc, Lauer & Postman 
(1994) found a dipole toward (I = 220°, b = -28°) 
±27°, inconsistent with the CMB dipole direction 
at 99.99% confidence and generated by mass con- 
centrations beyond lOOh^ 1 Mpc. An analysis of 
IRAS indicated a dipole direction within 13° of the 
CMB dipole that has not fallen to zero by 300/i _1 
Mpc (Rowan- Robinson et al. 2000). Kocevski et al. 
(2004) used an all-sky X-ray selected cluster cat- 
alog and found that the direction of the LG's pe- 
culiar velocity is well aligned with the CMB as 
early as the Great Attractor region 40/i _1 Mpc 
away. However, most of the Kocevski et al. (2004) 
dipole signal is attributed to the Shapley super- 
cluster 150/i _1 Mpc away, as well as a handful of 
massive clusters behind our galaxy. 

Using 56 SMAC clusters within 120/i _1 Mpc, Hud- 
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son et al. (2004) found a bulk flow of 687 ± 203 
km/s toward I = 260° ± 13°, b = 0° ± 11° that 
did not drop off with depth. They found that it 
could not be caused by the Great Attractor, but 
could be caused by the Shapley Concentration (at 
marginal significance). They argue that multiple 
data sets exclude convergence to the CMB frame 
by 60/1" 1 Mpc, and that at depths of 60 - 120/j" 1 
Mpc the flow is limited to 600 km/s. Using a 
new method of optimally weighting peculiar ve- 
locities, Watkins et al. (2009) compiled all major 
peculiar velocity surveys (including SMAC) and 
found them to be highly consistent. Within a re- 
gion of radius <~ 100ft. -1 Mpc they find a flow of 
407 ± 81 km/s toward I = 287° ± 9°, b = 8° ± 6° 
implying that 50% of the LG motion is generated 
beyond this depth. Extending this analysis Feld- 
man et al. (2010) find a consistent flow, 416 ± 78 
km/s toward I = 282° ± 11°, b = 6° ± 6°. These 
results are found to be in disagreement with the 
ACDM model with WMAP 5 year cosmological 
parameters at a high confidence level. 

Using the 2MASS Rcdshift Survey Lavaux et al. 
(2010) find that less than half of the amplitude of 
the CMB dipole is generated within ~40/i _1 Mpc, 
and that most of the amplitude of the dipole is re- 
covered by 120/i _1 Mpc, although the directions of 
the two flows do not agree. Colin et al. (2010) use 
Type la supernovae (SNe la) to find that there 
is a bulk flow of around 260 km/s at <~ 180/i _1 
Mpc, which disagrees with ACDM at the 1 — 2a 
level. However, at ~ 435/i _1 Mpc they find im- 
proved agreement between the SNe la data and 
the isotropic ACDM model. 

The kinetic Sunyaev-Zeldovich (SZ) effect was 
first used to place a limit on the bulk flow veloc- 
ity by Benson et al. (2003) using 10 clusters be- 
tween 300 - 600/1" 1 Mpc. They found no detec- 
tion but limited the flow in the direction of the 
CMB dipole to < 1410 km/s at 95% confidence. 
Kashlinsky & Atrio-Barandela (2000) found that 
by utilizing the kinetic SZ effect in WMAP data, 
flows as small as 200 km/s (or 30 km/s for the 
Planck experiment (ESA 2005)) could in principle 
be measured. Using the WMAP 5 year data Kash- 
linsky et al. (2008) (also Kashlinsky et al. 2009, 
hereafter KAKE) found a coherent dipole out to at 
least a distance of 300/i _1 Mpc. When all clusters 
in their sample beyond 300/i _1 Mpc are combined 
the dipole aligns well with the CMB dipole and is 
within 6° of the dipole found by Watkins et al. 
(2009). However, the magnitude of the KAKE 
dipole is considerably larger, estimated to be be- 
tween 600 — 1000 km/s and is detected on much 



larger scales with most of the signal coming from 
a region between 120 — 600/1^ Mpc away. Using a 
larger cluster sample Kashlinsky et al. (2010) find 
that the flow has constant velocity out to approx- 
imately 575/i -1 Mpc. However, using the same 
method as KAKE, Keisler (2009) do not find a 
significant detection of a bulk flow. They find that 
residual CMB in the filtered KAKE maps, that is 
correlated between the WMAP channels and not 
accounted for in the errors, decreases the signifi- 
cance of the cluster dipole. Atrio-Barandela et al. 
(2010) present an analysis of the errors in the clus- 
ter dipole measurement. They do not reproduce 
the errors in Keisler (2009) unless the monopolc 
and dipole are removed from the full sky CMB 
map instead of only in the region outside of the 
galactic mask, suggesting this as a reason for the 
null result of Keisler (2009). 

We measure the galaxy cluster peculiar velocity 
distribution using the WMAP 7 year data. The 
main difference between our method and that of 
KAKE and Keisler (2009) is the filter we use to 
suppress the CMB signal in the WMAP maps. We 
use a multi-frequency matched filter that optimizes 
detection of the cluster signal. As well as provid- 
ing a higher signal to noise measurement, our filter 
automatically accounts for CMB and noise correla- 
tions between the channels. Using simulations we 
find that if no attempt is made to reduce the ther- 
mal SZ signal in the maps, our results would be 
biased. We use a modified version of the matched 
filter that incorporates the thermal SZ spectrum, 
to reduce the thermal SZ signal by over an order of 
magnitude. In addition, we subtract the monopolc 
and dipole from the region outside of the galactic 
mask, and use an identical pipeline for the data 
and error analysis. 

We also measure the temperature monopole 
at the locations of galaxy clusters. The cluster 
monopole contains a contribution from the ther- 
mal and kinetic SZ effects. The kinetic SZ contri- 
bution is expected to be zero if there is no cluster 
monopole velocity, as we expect. In addition we 
measure the quadrupole, octupole, and £ = 4 and 5 
modes. Our method can be applied to data from 
the upcoming Planck experiment to achieve a re- 
sult with greater sensitivity (Mak et al. 2010). 

In section 2 we present all necessary theoretical 
background information. We describe the SZ effect 
and how it can be used to measure the bulk flow, 
and present the ACDM predictions for the velocity 
of the flow. In section 3 we describe the data that 
we use to make our measurement, and explain how 
we construct our cluster sample. In section 4 we 
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explain our method, describing how we model the 
SZ emission, how we construct filters to optimize 
cluster detection, and how we determine the am- 
plitude and direction of any bulk flow. In section 5 
we describe the level of contamination of thermal 
SZ, unresolved radio point sources, and galactic 
emission, which we determine using simulations. 
In section 6 we present our results and in section 7 
give our conclusions. 

Throughout this paper we assume a spatially 
flat ACDM cosmology with the WMAP 7 year cos- 
mological parameters (Jarosik et al. 2010). 

2. Theory 

2.1. SZ Clusters as Tracers of the Velocity 
Field 

The SZ effect can be used to measure the devia- 
tions of galaxy clusters from the Hubble flow (Sun- 
yaev & Zeldovich 1980). Approximately 1% of 
CMB photons traveling through galaxy clusters 
are scattered by electrons trapped in the gravita- 
tional potential of the cluster. If the galaxy cluster 
is moving with respect to the CMB rest frame, the 
scattered photons are red or blue shifted. This is 
the kinetic Sunyaev- Zeldovich (kSZ) effect (Sun- 
yaev & Zeldovich 1980) and results in a fractional 
temperature change in the radiation of 
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where dl is the line of sight distance through the 
cluster, v p is the cluster peculiar velocity, n e is the 
electron density and ot is the Thompson scatter- 
ing cross section. The quantity J n e ardl is the op- 
tical depth to Thompson scattering, r, with typical 
value ~ 0.01 for the clusters we observe. By observ- 
ing the temperature increment or decrement in the 
direction of a galaxy cluster, the line of sight clus- 
ter peculiar velocity can be estimated. For clus- 
ters in our sample we expect ATksz/T ~ 10 -5 . 
To date no detection of a cluster velocity has been 
made using the kSZ effect. 

A bulk motion of mass in the universe will cause 
a dipole pattern in temperature at the cluster po- 
sitions. The kSZ signal from each cluster that is 
part of a bulk motion is given by, 
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where r is the optical depth of the cluster and 6 
is the angle between the cluster position and the 
bulk flow direction. Cluster motions with more 



structure than monopole and dipole terms can be 
described by decomposing the cluster velocity field 
into spherical harmonics. 

The electrons within massive clusters have tem- 
peratures that are typically a few keV. The ther- 
mal motion of these electrons causes an increase 
in the energy of the scattered photons, which is 
the thermal Sunyaev-Zeldovich (tSZ) effect (Sun- 
yaev & Zeldovich 1970). The fractional tempera- 
ture change in the radiation in the non-relativistic 
limit is 



77 k BT e 

dl ^n e aT (3) 

m e cr 



where x = hv jUsTcus, v is the radiation fre- 
quency, dl is the line of sight distance through the 
cluster and T e is the electron temperature. To de- 
termine the importance of relativistic corrections 
to equation 3, we have calculated the thermal SZ 
spectrum of all 736 clusters in our cluster sample 
(we describe the cluster modeling we use in sec- 
tion 4.2), using equation 3 to calculate the tSZ 
emission and also including relativistic corrections 
using the model of (Nozawa et al. 2000). We com- 
bine the spectra, weighting the spectrum from each 
cluster with the same weight that we use when we 
measure the bulk flow signal (we describe how we 
calculate the bulk flow signal in section 4.4). We 
find that the difference between the flux calculated 
using equation 3 and the flux from the relativistic 
formula is less than 5% at all of the frequencies we 
use, and so for our purposes we can safely ignore 
the relativistic corrections. 

For typical clusters AT TSZ /T ~ 10~ 4 at 90 
GHz, and the ratio of the kinetic signal to the ther- 
mal signal is (assuming the clusters are isother- 
mal), 
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The thermal signal is therefore a significant con- 
taminant to any measurement of the kinetic signal. 
However, since there are not expected to be intrin- 
sic large scale moments in the tSZ signal, the tSZ 
dipole amplitude is expected to be smaller than 
the monopole amplitude. 

Unlike the tSZ signal, the kSZ signal has an 
identical frequency spectrum to the CMB. The 
kSZ signal can therefore only be distinguished from 
CMB fluctuations by the different spatial proper- 
ties of the two signals. The CMB power peaks on 
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degree scales, whereas the SZ emission from galaxy 
clusters typically varies on arcminute scales. The 
signal we detect is smeared out by the instrumen- 
tal beams of the experiment, 30'. 6 in the WMAP 
Q band, 21' in the V band and 13'. 2 in the W 
band. 

2.2. Expected Signal For a Generic Veloc- 
ity Tracer 

On scales larger than <~ 10/i _1 Mpc, fluctuations 
in the matter density of the universe are Gaussian 
distributed (Rimes & Hamilton 2005, and refer- 
ences therein). The inhomogeneities in the matter 
density cause galaxy clusters to have peculiar ve- 
locities, with typical values of ~ 300 km/s (e.g., 
Bhattacharya & Kosowsky 2008). If a sample of 
galaxy clusters within a region is chosen, that sam- 
ple will have a non-zero bulk motion due to the 
non-uniform matter density on even larger scales. 
Although the direction of this bulk motion is not 
determined by the ACDM model, the rms ampli- 
tude can be calculated given a set of cosmological 
parameters. We now calculate this amplitude. 

We can estimate the expected cluster bulk flow 
velocity by expanding the line of sight peculiar 
velocity distribution in spherical harmonics. The 
rms bulk flow velocity will be the power in the 
first multipole. We follow the derivation for the 
density distribution given in Peebles (1973). We 
decompose the line of sight peculiar velocity into 
spherical harmonics, Yg m (r), with the amplitude 
of the £,m mode given by: 
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where r is the comoving radial distance, </> is the 
comoving number density of objects in the sample, 
v is the object peculiar velocity and fi r is the solid 
angle. We approximate <f)(r) by an isotropic func- 
tion, which we estimate by calculating the number 
density of clusters in our sample within radius r. 
We normalize 6 such that: 
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where R is the comoving distance within which 
the ai m are calculated. Figure 1 shows the redshift 
distribution of our cluster sample used to calculate 

<t>. 

The power in multipole £ is defined as Cg = 
(\de m \ 2 ) where the average is over all values of m 



within multipole I. If the line of sight peculiar ve- 
locity over the whole sky is a Gaussian random 
field, as expected in the ACDM model, then it 
is completely described by C'e. The power in the 
dipolc {£ = 1) is (the details of the derivation arc 
given in appendix A): 
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where P{k) is the matter power spectrum, / = 
(a/D) dD/da, a is the scale factor, D is the growth 
function, Hq is the Hubble constant and j\ and 
ji are spherical Bessel functions. Since the clus- 
ters we observe are all at redshift less than one, 
we approximate / by its value today, taking it to 
be equal to £l m 6 - We obtain a prediction of the 
rms dipole velocity shown in figure 2. Since our 
cluster sample contains few clusters with redshifts 
greater than ~ 0.3, the bulk velocity of our en- 
tire cluster sample is largely determined by clus- 
ters with redshifts less than this. The effect of 
the selection function is therefore to increase the 
expected dipole velocity in shells extending to red- 
shifts greater than 0.3. The shaded area in figure 2 
is the uncertainty from cosmic variance. The ad- 
ditional uncertainty on C\ from sample variance 
is inversely proportional to the number of clusters 
that are used to measure the dipole. For redshift 
shells extending beyond z = 0.05 our cluster sam- 
ple contains over one hundred clusters and so the 
sample variance in figure 2 will be more than ten 
times smaller than the cosmic variance. We there- 
fore do not include it in the figure. 

A measurement of a kSZ bulk flow with an am- 
plitude larger than that expected from cosmic vari- 
ance would be an important result if confirmed, re- 
quiring modifications of inflation to explain it (e.g., 
Turner 1991; Kashlinsky et al. 1994). A measure- 
ment of a kSZ monopole with a velocity greater 
than that expected by sample variance would be a 
violation of the Copernican principle, which states 
that we do not live in a specially favored place in 
the universe. However, a measurement of the kSZ 
monopole is more susceptible to contamination. 
The thermal SZ signal, as well as the radio point 
source signal from clusters, are both expected to be 
isotropic, with any dipole signal coming from sam- 
ple variance. Any dipole signal from tSZ or radio 
point sources is therefore expected to be weaker 
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Fig. 1. — Redshift distribution for the cluster sam- 
ple used in this paper. 
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Fig. 2. — Expected bulk flow velocity in ACDM 
cosmology with a selection function calculated 
from our cluster sample in redshift shells extend- 
ing from z = to the specified x-axis value. The 
dashed line shows the result with a uniform se- 
lection function. The dark and light shaded areas 
are the 68% and 95% confidence limits from cosmic 
variance. 



than any corresponding monopole signal. The ra- 
dio source monopole at the cluster locations could 
be removed by observing and subtracting the emis- 
sion, or by utilizing the different spectra of the ra- 
dio source and kSZ signals. For thermal SZ the 
monopole can be suppressed spectrally. 

3. Data 

We have searched for a cluster dipole in the 
WMAP 7 year maps 1 (Jarosik ct al. 2010). We use 
the foreground reduced maps produced by WMAP. 
These were generated by removing a foreground 
model from each 'unreduced' map (Hinshaw et al. 
2007). Both the K and Ka band maps were used 
to produce the galactic foreground model, and so 
there are no K or Ka band foreground reduced 
maps. The maps are pixelized in the Healpix for- 
mat (Gorski et al. 2005) with a pixel size of 7'. 
The CMB dipole has been removed from each map. 
Removing the CMB dipole from the maps does 
not significantly change the cluster dipole that we 
search for. We have tested the effect of remov- 
ing the CMB dipole by simulating kSZ maps (our 
simulation procedure is described in section 4.2), 
adding a CMB dipole signal, and then removing 
it using a least squares fit to all map pixels. We 
find that the amplitude of the cluster dipole in 
the resulting Q band maps is changed by less than 
1%, and the direction of the dipole is changed by 
less than 25'. In the V and W band channels the 
change is smaller due to a stronger cluster signal 
relative to the CMB dipole signal. 

We use the two Q (41 GHz), two V (61 GHz) 
and four W (94 GHz) band maps from each dif- 
ferencing assembly. We exclude the K and Ka 
band maps from the analysis since no foreground 
reduced maps are available. The galactic syn- 
chrotron and extra-galactic radio emission has a 
spectral index a > (with source flux, S <~ v~ a ), 
and the emission is stronger in the K and Ka band 
than in the Q, V or W bands. Although the instru- 
ment noise is lower in the K and Ka bands than in 
the Q, V or W bands, the cluster signal is weaker, 
due to dilution by the larger beam FWHM. In- 
cluding the K and Ka band channels therefore only 
increases the signal to noise of the cluster dipole 
measurement by ~ 10% and so we exclude them 
from our analysis. 

The noise per pixel is ~ 65^K (Q band), ~ 80^K 
(V band), and - 137^K (W band). The noise is 



The WMAP maps, beam functions and galactic masks are 
available at http://lambda.gsfc.nasa.gov/product/map/ 
current/m_products . cf m 
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lower near the ecliptic poles due to the WMAP 
scan strategy. The beam FWHM are 30'. 6 (Q 
band), 21' (V band), and 13'. 2 (W band). We 
generate noise simulations using the prescription 
described in Limon et al. (2006). The noise is as- 
sumed to be uncorrelated between pixels, which is 
a good approximation for the temperature maps 
that we use (Limon ct al. 2006). We use the 
WMAP beam transfer functions, which are the 
square roots of the window functions used for 
the power spectrum analysis, to construct our 
matched filters and smooth our simulated maps. 
The WMAP extended temperature analysis mask 
is used to exclude the galaxy and other known 
regions of high foreground emission. Any cluster 
that lies behind this mask is not included in our 
cluster sample. 

3.1. The X-ray Selected Cluster Sample 

Following KAKE we derive our cluster sam- 
ple from the REFLEX (Bohringer et al. 2004), 
BCS (Ebeling et al. 2000b), Extended BCS (Ebel- 
ing et al. 2000a) and CIZA (Ebeling et al. 2002, 
Kocevski et al. 2007) cluster catalogs. These cat- 
alogs contain 447, 206, 107, and 130 clusters re- 
spectively. After removing 20 overlapping clusters 
we obtain a sample of 870 clusters. After remov- 
ing clusters whose center lies behind the WMAP 
galactic mask we find 736 clusters (when using the 
WMAP KpO galactic mask that was used in the 
KAKE analysis we find 771 clusters in the sam- 
ple). The redshift distribution of the cluster sam- 
ple is shown in figure 1. We convert the luminosi- 
ties in the BCS, BCSe and CIZA catalogs to the 
values for a ACDM universe. For the REFLEX 
sample we use a catalog with the luminosities al- 
ready converted to an h = 0.7, ACDM cosmology 
so no correction is necessary for these clusters. 

To derive average electron density and temper- 
ature profiles for our cluster sample we use the 
Archive of Chandra Cluster Entropy Profile Ta- 
bles (ACCEPT, Donahue et al. 2006; Cavagnolo 
ct al. 2009). This is a sample of 239 clusters ob- 
served by the Chandra X-ray Observatory (Wcis- 
skopf et al. 2000), with electron density and tem- 
perature profiles for each cluster derived from a 
deprojection analysis. 145 of the ACCEPT clus- 
ters overlap with our sample. 

3.2. External Simulations 

To test for contamination from galactic emis- 
sion, thermal SZ and unresolved radio sources, we 
use simulated maps produced using the Planck Sky 



Model 2 (PSM). 

The PSM is a set of programs and data used 
for the simulation of full sky microwave maps that 
includes galactic synchrotron, dust and free-free 
emission, kinetic and thermal SZ, radio and in- 
frared point sources and CMB. The PSM simu- 
lations do not include a kinetic SZ monopole or 
cluster bulk flow velocity. We use one realization 
of the sky generated at the WMAP Q, V and W 
band center frequencies. 

To test the effect of radio point sources we use 
100 simulated maps of the unresolved radio back- 
ground produced by Colombo & Pierpaoli (2010). 
They combine sources from the NRAO-VLA Sky 
Survey (NVSS) catalog (Condon et al. 1998) with 
higher frequency surveys to extrapolate the radio 
point source emission at 1.4 GHz to the WMAP 
frequencies. The radio source flux distribution in 
the maps at each frequency therefore reflects the 
probability distribution for the flux of each NVSS 
source. 

4. Method 

4.1. Outline of the Method 

The bulk motion of many galaxy clusters rela- 
tive to the CMB rest frame creates a signal in the 
WMAP temperature maps that has a dipole pat- 
tern at the location of clusters. We now describe 
the method we use to recover the amplitude and 
direction of any dipole motion using the WMAP 
CMB maps. 

In addition to kSZ the dipole signal at the posi- 
tion of clusters includes contributions from CMB, 
instrumental noise, tSZ, galactic emission and in- 
frared and radio point source emission. The first 
step is to filter the maps to enhance the kSZ sig- 
nal relative to the other terms. We use two filters, 
one of which enhances the cluster signal relative 
to the CMB and instrument noise, the other filter 
also removes the thermal SZ term. The galactic 
term is suppressed by masking the WMAP maps 
with the WMAP extended temperature analysis 
mask, and only fitting for a cluster dipole outside 
of this region. In section 5.2 we give upper lim- 
its on the radio point source contamination and 
in section 5.3 we give upper limits on the galac- 
tic signal. Using the PSM simulations we find a 
negligible contribution from infrared point source 
emission. 



2 http : //www . ape . univ-paris7 . f r/APC_CS/Recherche/ 
Adamis/PSM/psky- en . php 
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The next step is to calculate the monopole and 
dipole at the locations of the ROSAT observed 
clusters in the filtered maps. We use simulated 
kinetic SZ maps to convert the amplitude of the 
dipole at the cluster locations to a bulk flow veloc- 
ity using a method we describe in section 4.5. Our 
simulated kinetic and thermal SZ maps include sig- 
nal only from the ROSAT observed clusters that 
we use, and will be described in section 4.2. 

There are two ways we could calculate the clus- 
ter velocity dipole. We could calculate the opti- 
cal depth through each cluster and use equation 1 
to estimate the velocity of each cluster. We could 
then fit for a dipole in the cluster velocity distribu- 
tion. Alternatively, we could calculate the dipole 
at the cluster positions in the WMAP maps and 
convert the dipole signal to a velocity using an 
estimate of the average optical depth of the clus- 
ters in our sample. The two methods differ in the 
weight given to each cluster in the dipole fit. We 
expect the signal to noise of the latter method to 
be higher since all clusters are given equal weight 
in the fit, whereas in the former method clusters 
with a larger optical depth (and hence a greater 
kSZ signal in the map) are given less weight. We 
follow KAKE and choose the latter method. We 
split the cluster sample into different redshift bins 
(some of which overlap) and calculate the cluster 
dipole in each. 

To calculate the errors on the monopole, dipole 
and higher order modes, we generate 100 realiza- 
tions of the CMB, convolve the maps with the 
beams from each of the eight channels and add 
instrument noise using the procedure described 
in Limon ct al. (2006). We pass these maps 
through our pipeline to find the distribution of 
monopole, dipole and higher order modes. 

4.2. Cluster modeling 

Our measurement of the WMAP temperature 
dipole does not depend on any cluster simulations. 
However, in order to convert the cluster tempera- 
ture dipole measurement into a velocity dipole we 
require simulated maps of the kSZ emission from 
our cluster sample. In addition, we require tSZ 
realizations to estimate the level of tSZ contami- 
nation in our results. 

The clusters in our sample have typical angular 
sizes of ~ 1', which is small compared with the 
WMAP beam FWHM, 30'.6 in Q band and 13'. 2 
in W band. For this reason we model the clusters 
as point sources convolved by the WMAP beams. 



4-2.1. Optical Depth Determination 

In order to generate kSZ realizations with a 
simulated bulk flow velocity we require the opti- 
cal depth to Thompson scattering for every clus- 
ter in our sample. We calculate this using aver- 
age electron density and temperature profiles that 
we derive from the ACCEPT catalog (described in 
section 3.1). 

To calculate the electron density and tempera- 
ture profiles we select 145 ACCEPT clusters that 
overlap with our sample. We average the electron 
density and temperature profiles of the overlap- 
ping clusters within radial bins. We use bins with 
a smaller size nearer to the center of the cluster 
where the error on the electron density and tem- 
perature is lower (the error is provided in the cat- 
alog). The mean bin size is 10 kpc, and the largest 
radial bin extends to 1 Mpc from the cluster cen- 
ter. Each profile is normalized so that it has unit 
value at the center of the cluster. The densities 
and temperatures of all of the clusters are then 
averaged within each radial bin. We compute the 
kSZ and tSZ signals using these profiles, giving the 
profiles a different normalization for each cluster, 
as described in the following paragraph. Extrap- 
olating the density profile to larger radii using a 
/3— model (Cavaliere & Fusco-Femiano 1976), we 
find a negligible contribution to the optical depth 
from beyond this region. The profiles are not sig- 
nificantly changed when all of the clusters in the 
ACCEPT sample are used to calculate them. 

Since we use the same electron density profile 
for each cluster, we only need to calculate the nor- 
malization to that profile to estimate the optical 
depth. The optical depth is calculated by integrat- 
ing the electron density along the line of sight, 

where R = cLaO is the distance from the cluster 
center perpendicular to the line of sight, oIa is the 
angular diameter distance, ctt is the Thompson 
scattering cross section, f e (r) is the normalized 
electron density profile and N e is the normaliza- 
tion to that profile. 

We calculate the normalization to the elec- 
tron density profile by requiring that the cluster 
bolometric luminosity in the ROSAT catalogs is 
equal to the luminosity from Brcmsstrahlung emis- 
sion, which depends on N e (Gronenschild & Mcwe 
1978), 
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4-2.2. kSZ and tSZ Simulations 



32tt {2irk B y Q 6 ^ 

L =-5- tj-Z gN e Ni 

3 V 3m e J m e c A h ^ 

J d 2 rAnr 2 f e (r)fi(r)T(r^ erg/s 

where ks is the Boltzmann constant, m e is the 
electron mass, Q is the electron charge, c is the 
speed of light, h is Planck's constant, Z is the num- 
ber of protons per nucleus, g is the frequency av- 
eraged Gaunt factor (of order unity), T{r) is the 
electron temperature, f e {r) and fi(r) are the nor- 
malized electron and ion density profiles (which we 
take to be equal) and N e and JV, are the normaliza- 
tions to those profiles. We normalize the tempera- 
ture profile such that the X-ray emission weighted 
temperature is equal to T\, which we calculate 
from the X-ray luminosity of each cluster using 
the relation T x = (2.76 ± 0.08)L^ 33±001 (White 
et al. 1997). 

Using this procedure we find that the mean cen- 
tral optical depth of our cluster sample is (4.9 ± 
0.9) x 10~ 3 . The mean central optical depth of 
the 145 clusters that overlap with the ACCEPT 
sample is t = (6.1 ± 1.0) x 10~ 3 which is in agree- 
ment with the value calculated using each individ- 
ual electron density profile in the ACCEPT cata- 
log, r = (6.7 ± 0.9) x 10~ 3 . This gives us con- 
fidence that we can calculate accurate values of 
the optical depth for clusters that are not in the 
ACCEPT sample. The dominant source of uncer- 
tainty in our optical depth calculation is the cluster 
X-ray luminosity, from the X-ray catalogs. Propa- 
gating all errors through to the optical depth gives 
uncertainties of ~ 15%. These errors are propa- 
gated through to the simulated kSZ maps which 
allows us to calculate the precision of our method. 

We have not accounted for any anisotropy in 
the optical depth distribution of our clusters or 
accounted for the different methods used to calcu- 
late the X-ray luminosities in the different X-ray 
catalogs we use. Both of these steps would require 
an independent analysis of the X-ray data. How- 
ever, since our dipole fit does not depend on the 
X-ray luminosities our null detection is unaffected. 
Any error in the X-ray luminosities will only affect 
the conversion of the dipole amplitude from /xK to 
km/s. Since the X-ray luminosities are used to 
calculate the optical depth any asymmetry in the 
optical depths will be accounted for when convert- 
ing from fiK to km/s. 



The kSZ signal is calculated from equation 1. 
We give each cluster in our simulated maps a bulk 
velocity using equation 2 and a random velocity 
drawn from a Gaussian distribution, with vari- 
ance (e.g., Gorski 1988) 

<r 2 v = ^ J dk P(k) \W(kR)\ 2 (10) 

where |VF(fci?)| 2 is the Fourier transform of the 
window function, and again we use / ~ O^ 6 . To 
select cluster scales we use a top hat window func- 
tion with R = Wh- 1 Mpc, and find a v = 272 
km/s. We ignore correlations between the ve- 
locities of different clusters, which we expect to 
be small for our cluster sample (eg. Gorski 1988; 
Bhattacharya & Kosowsky 2008). 

The expected tSZ signal is calculated by inte- 
grating equation 3 through the cluster using the 
same electron density and temperature profiles cal- 
culated in section 4.2.1. Since WMAP is sensitive 
enough to detect the tSZ signal in stacked images 
we have checked that the average tSZ signal in 
our simulations agrees with WMAP observations. 
Following the method in Atrio-Barandela et al. 
(2008) (also Diego & Partridge 2010), we have 
stacked maps centered on each cluster to produce 
a coadded map in the Q, V and W bands from the 
WMAP maps and our simulated tSZ maps. We 
find that the residual after subtracting the two is 
consistent with CMB and instrument noise. 

In table 1 we show the simulated thermal SZ 
monopole in maps filtered by the matched filter 
we use in our analysis. We compare this to the 
monopole found in the WMAP maps filtered by 
the same filter. Although the monopole in our 
simulated maps is biased high it is within 2a of 
the measured monopole in all except the two low- 
est redshift shells. This gives us confidence that 
our optical depth calculation is reasonable. 

4.3. Filters 

The distribution of our cluster sample on the 
sky is not isotropic and so measurements of the 
peculiar velocity multipoles will be contaminated 
by higher order CMB multipoles. The higher order 
CMB multipoles couple to the signal we are trying 
to measure through the mask we use to cover non- 
cluster regions of the WMAP maps. To reduce the 
contamination from the CMB we filter the maps. 

Two filters are used, and described here. The 
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Table 1 

Simulated thermal SZ monopole and measured monopole. 



Zmin z m ax N c \ Simulated tSZ Monopole, ao [/iK] Measured Monopole, ao [/iK] 



0.0 


0.04 


95 


-78.7 


-22 ± 33 


0.0 


0.05 


139 


-84.2 


-41 ± 27 


0.0 


0.06 


192 


-97.0 


-62 ± 24 


0.0 


0.08 


294 


-97.9 


-63 ± 18 


0.0 


0.12 


445 


-99.5 


-73 ± 13 


0.0 


0.16 


546 


-103.4 


-81 ± 12 


0.0 


0.20 


619 


-109.1 


-85 ± 12 


0.0 


1.0 


736 


-123.9 


-91 ± 11 


0.05 


0.30 


578 


-128.6 


-99 ± 12 


0.12 


0.30 


271 


-152.2 


-112 ± 17 



first is an optimal matched filter that removes low- 
order CMB multipoles by utilizing information at 
locations other than the cluster position to distin- 
guish a CMB multipole from a kSZ multipole. The 
second filter uses both spatial information and the 
spectral difference between the thermal and kine- 
matic SZ effects to remove the contaminating tSZ 
signal. 

An optimal filter is constructed for each channel 
such that when the filtered maps from each channel 
are combined, the kSZ signal to noise ratio is max- 
imized. In order to construct our filters we require 
a matrix of cross spectra, C e = C^ MB + C£ oisc , 
that describes the statistical properties of the maps 
to be filtered. We assume that the noise in map 
v, n v (0), is a homogeneous and isotropic random 
field. The cross power spectrum of maps v\ and 

v 2 is equal to (rii mtI/1 n} m ^ , where n im ^ arc the 
coefficients in the spherical harmonic expansion of 
n v {6), and the average is over to. We generate Ce 
by averaging the cross power spectra of 100 CMB 
and noise realizations. 

The maps in each channel are filtered and com- 
bined to produce the spherical harmonic coeffi- 
cients of the filtered map as follows, 

"channels 

hm= ®e a e™ (11) 

i=l 

where a\ m are the spherical harmonic coefficients 
of the WMAP channel i map and $\ is the filter 
for channel i. We include the factor of y/Air/21 + 1 
from spherical convolution in The first filter 
we use has a filter function of the form (McEwen 
et al. 2008): 



r 1 

= -^-B £ (12) 

where m indicates that the filter is a matched fil- 
ter, 7 = (1/ripix) J2e(^ + l)Bj , C7 1 B £ , n pix is the 
number of map pixels and (47r/np ix )B| is an es- 
timate of the kSZ power spectrum from a single 
cluster, in each channel (we assume the beam con- 
volved sources are symmetric). Since the WMAP 
beams are much larger than galaxy clusters we 
approximate by the WMAP beam function. 
We have calculated the filter using the WMAP 
beam profile convolved with an average cluster 
profile and also a cluster profile modeled as a 
f3— model (Cavaliere & Fusco-Femiano 1976) and 
find no significant difference in the results. 

A simple modification to this filter allows the 
thermal SZ bias to be removed (e.g., Herranz et al. 
2005; Schafer et al. 2006, for the flat and curved 
sky cases respectively), 

$ " = %" {aBe ~ (13) 

where u indicates that the filter is an unbiased 
filter, Fg is an estimate of the thermal SZ sig- 
nal in each channel, A = 017 — j3 2 is a normal- 
ization factor, a = (l/n pix ) J2e(2£ + l)FjC^ 1 F^, 
= (1/npix) Ee( 2i + lJBjC^F* and again 7 = 
(l/ n pix) l)Bj C^Bf . We are using maps 

with units of thermodynamic temperature and so 
F t = [x(e x + l)/(e x - l)-4.0] B e . 

Figures 3 and 4 show the functions $™ and 
for each channel. These filters have very different 
behaviors and we now discuss why this arises. 

The matched filter (hereafter MF) suppresses 
the input map at multipoles less than <~ 300 (large 
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Fig. 3. — Multi- frequency matched niters for the 
WMAP channels. 
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Fig. 4. — Unbiased multi-frequency matched fil- 
ters. A discussion of the filter shape is given in 
the text. 



scales) due to the strong CMB signal at these mul- 
tipoles. The filter function also suppresses the map 
at multipoles greater than ~ 1000 (small scales), 
where the expected kSZ signal to noise ratio is 
low. The W band filter is non-zero up to the high- 
est multipoles since the W band beams are smaller, 
and so the signal extends to higher multipoles. Be- 
tween multipoles 300-1000 the signal to noise ratio 
is highest and so these scales are retained in the 
filtered map. The bumps near multipoles 500 and 
700 in the MF correspond to the troughs in the 
CMB spectrum— the filter amplifies scales where 
the contaminating CMB emission is weaker. 

The tSZ bias removing filter (hereafter UF) in 
figure 4 is less intuitive. At low multipoles the 
CMB signal in each channel is almost identical be- 
cause each channel observes the same CMB sky, 
and at low multipoles the beam functions are all 
close to unity. At low multipoles the CMB can 
therefore be removed by subtracting the Q and W 
band maps. Since there are two Q band channels 
and four W band channels the CMB is removed by 
giving the Q band maps double the weight of the 
W band maps in the subtraction. This explains 
why the two Q band filters have an absolute value 
double that of the four W band filters at low mul- 
tipoles. This is reflected in equation 13 by the 
off-diagonal elements being almost equal to the di- 
agonal elements in the cross-power spectra matrix 
at low multipoles. If the CMB sky in each chan- 
nel were different then the filter would look similar 
to the MF but with different amplitudes for each 
channel. In figure 5 we show what the tSZ bias 
removing filters would look like for an experiment 
with four frequency channels: 100 GHz, 143 GHz, 
217 GHz and 353 GHz, with beam FWHM of 10', 
7', 5' and 5' respectively, and white noise levels: 
25, 15, 25 and 75 ^K/K-arcmin. The fact that 
there are no peaks in the WMAP filters is a reflec- 
tion of the limited frequency coverage and large 
beams. 

Since the UF combines channels in a way that 
removes the tSZ signal, the filters do not increase 
the kSZ signal to noise ratio as much as the MF. 
The signal to noise of the cluster temperature 
dipole measurement is therefore lower. 

Figure 6 shows the spectra of a simulated map 
containing only the kSZ signal at the ROSAT clus- 
ter locations with and without a bulk flow compo- 
nent of 1000 km/s, and convolved by the WMAP 
beams. The maps with and without a bulk flow 
both have spectra that are identical to the WMAP 
beam function at high multipoles. The map with 
the bulk flow component contains more power at 
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all multipoles as well as additional power at the 
dipole and octupole modes. Although the dif- 
ference between the spectra with and without a 
bulk flow looks predominantly like an amplitude 
shift, the bulk flow signal is recoverable by in- 
formation contained in the spherical harmonic m 
modes, which are averaged to create the spectra. 
Our filters are designed to optimize the detection 
of sources shaped like the WMAP beams. The 
similarity of the simulated kSZ spectrum with a 
bulk flow component to the beam profile indicates 
that our filter will have the desired effect of in- 
creasing the signal to noise of the cluster dipole 
measurement. 

Figures 7 and 8 show the WMAP maps after 
filtering with the MF and UF respectively. Both 
filters have suppressed the CMB signal and the 
only visible structure in the map is the low noise 
region around the ecliptic poles. The maps are 
noise dominated; the larger values of the map in 
figure 8 reflects the lower signal to noise ratio ob- 
tained with the UF. 

Since we are only calculating the dipole in a 
small region at the center of the clusters, we could 
filter the maps only in this region. This would 
save computation time but care would need to be 
taken to avoid ringing effects around the edge of 
the mask. For this reason we filter the maps in all 
regions outside the WMAP galactic mask. 

4- 3.1. Wiener filters 

We now compare our filter to the Wiener filter 
used by KAKE. A Wiener filter (hereafter WF) 
minimizes the squared difference between the fil- 
tered map and an estimate of the signal. Although 
a multi-frequency WF can be constructed, we fol- 
low KAKE and create a filter for each map that 
is designed to suppress the CMB component. The 
KAKE filter is (Kashlinsky et al. 2009): 

$ w = Ce ll*yz±2 — £i ( M ) 

W / /sky 

where C^ MB is the CMB power spectrum, Be is 
the WMAP beam function, C| ky is the spectrum 
of the WMAP map estimated outside of the galac- 
tic mask and / s k y is the fraction of the sky out- 
side of the mask. The full sky spectra could alter- 
natively be calculated by deconvolving the mask: 

C unmasked l\ /f — 1 /^masked /^masked / £ 

I — lui' 1V1 W W ~ W //sky, 

where Mu> is the multipole mixing matrix which 
accounts for the cut sky and is calculated from 
the galactic mask (see eg. appendix A of Hivon 
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; 

Fig. 5. — Thermal SZ bias removing filters for an 
experiment with four frequencies: 100 GHz, 143 
GHz, 217 GHz and 353 GHz, beam FWHM of 10', 
7', 5' and 5', and white noise levels 25, 15, 25 and 
75 jtxK/K— arcmin. 




Fig. 6. — Spectra of a simulated kSZ map (with 
signal only at the locations of our cluster sample) 
convolved with the WMAP beams. The simula- 
tions were performed both with (top 3 spectra) 
and without (bottom 3 spectra) a bulk flow com- 
ponent of amplitude 1000 km/s. The dashed lines 
are the WMAP beam functions, scaled to the sig- 
nal amplitude. 



11 




Fig. 7. — Sum of the eight WMAP maps in the Q, 
V and W bands filtered by the matched filters in 
figure 3. The map is noise dominated. 




ct al. (2002) for details). The KAKE filter we use 
is shown in figure 9. We have determined that 
this is the same Wiener filter that was used in the 
KAKE analysis and give details of our pipeline in 
appendix B. The KAKE filter is constructed from 
the spectrum of the map that is to be filtered, and 
maps filtered with it suffer from low multipole fluc- 
tuations, caused by cosmic variance. We create our 
Wiener filter in such a way that it does not suffer 
from this problem (Tegmark & Efstathiou 1996), 



C noise 
I 



//sky 



rtCMB r>2 i /-^noisc / f 
W n l + 1 /sky 



(15) 



where C" OISC is the power spectrum of the instru- 
ment noise plus all foreground components in the 
masked map. This filter is not normalized and 
so pixels in the filtered map are not equal to the 
estimated amplitudes of the kSZ signal in the un- 
filtered map pixels. Since we model our clusters 
as point sources convolved by the beam the nor- 
malization factor would be n v \^j J2e(^ + ^)Bi$g- 
The normalization is instead applied later when 
the dipole amplitude is converted from Kelvin to 
km/s. Our Wiener filter is shown in figure 10. 

Since the noise in the WMAP maps has a flat 
spectrum and the CMB signal decreases at high 
multipoles, the WF acts like a high pass filter for 
each of the WMAP channels. At low multipoles 
the shape of this filter is similar to our MF (al- 
though the relative normalizations of the channels 
is not), but at high multipoles the shapes differ 
since the WF does not take into account the fact 
that the cluster signal decays exponentially with 
increasing multipole, while the noise spectrum re- 
mains flat. We find that both our Wiener filter 
and the KAKE filter suppress the cluster signal 




10.3 mK 



Fig. 8. — Sum of the eight WMAP maps in the 
Q, V and W bands filtered by the tSZ removing 
filters in figure 4. The values in this map are larger 
than in figure 7 since the map is noise dominated 
and the unbiased filters give a lower signal to noise 
measurement of the kSZ signal. 




Fig. 9.— Wiener filter used by KAKE. 
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approximately equally as shown in figure 11. At 
the multipoles where the signal peaks (~ 400 in 
the V band) the two filters suppress power by a 
similar amount. We find that our Wiener filter 
reduces the noise in the a\ m by a factor of ~ 2 rel- 
ative to the KAKE filter by reducing CMB noise 
at multipoles less than ~ 100 that is correlated 
between the WMAP channels (we give a more de- 
tailed comparison of the filters in appendix B) , and 
is therefore more sensitive to the cluster signal. 

For the matched filter and unbiased filter we cal- 
culate the cluster dipole in the pixels that overlap 
with the centers of each cluster. For our Wiener 
filter we instead use a 15' aperture around each 
cluster, in order to increase the signal to noise of 
the measurement. We have tried fitting for the 
dipole using different aperture sizes out to 30' as 
well as different weighting schemes for the pixels 
within the aperture: uniform weighting, weight- 
ing by the inverse noise variance and weighting by 
the expected cluster signal to noise ratio. The re- 
sults presented in table 4 have a uniform weight- 
ing within the aperture, which was chosen to give 
a good signal to noise measurement while keeping 
our results free from any bias caused by incorrect 
signal or noise estimates. 

The WF that we use is not a multi-frequency 
filter (although in principle a multi-frequency WF 
could be used) and so there is a filtered map for 
each of the WMAP channels. We filter each map 
and compute the cluster dipole in each. We then 
average the dipoles from the different channels 
with uniform weights: 



01, 



1 8 



(16) 



The signal to noise could be increased by 
weighting the results from each channel by the 
expected signal to noise ratio of the kSZ signal. In 
this scheme the weights given to each channel are 
approximately equal— the noise at higher frequen- 
cies is compensated for by the decreased beam 
FWHM, and therefore larger cluster signal. 

4.4. Dipole Fitting Procedure 

To measure the cluster dipole we perform a 
weighted least squares fit to a dipole function at 
the locations of the clusters in the filtered map. 
We also fit for the monopole and, separately, we 
perform the fit for modes up to I = 2, and up to 
£ = 3,4,5. Our dipole fitting method is based on 
the Healpix IDL procedure remove-dipole (Gorski 
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Fig. 10.— Wiener filter for each of the WMAP 
channels. 
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Fig. 11. — Suppression of the kSZ signal by the 
beam and filter in the V band. 
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et al. 2005). We fit for a vector of coefficients, 
(3, of the real spherical harmonics using the least 
squares formula, 

p = (X T WX)- 1 X T Wy (17) 

where y is the data map, W is a diagonal matrix 
with diagonal values equal to a weights map, and 
X is a matrix giving the contribution of the fitting 
function to each pixel. If a monopole and three 
dipole coefficients are fit to the data, then X is 
an n p ix x 4 matrix. We choose the weights map 
to have non-zero values only in the central cluster 
pixels, and equal to the inverse noise variance of 
the filtered map. The noise variance is calculated 
from 100 filtered WMAP CMB and noise realiza- 
tions. The vector (3 are the coefficients of the real 
spherical harmonics, Ri m , defined by 

r Y m if ™ = o 

Rim = \ ^ (Y em + (-l) m Y e _ m ) if to > 
{^(Ye- m -(-l) m Y em ) ifTO<0 

(18) 

The matrix X T WX is a mixing matrix that 
couples different spherical harmonic modes to- 
gether. The effect of the mask W is that a least 
squares fitting process can then be used on the 
masked map to determine the dipole from clusters 
alone. Caution is needed because information is 
lost when the map is masked. Consequently if too 
few modes are fitted, power will 'leak' from higher 
order modes into the fit parameters, corrupting the 
result. 

Since the kSZ signal is proportional to the op- 
tical depth, in principle the signal to noise ratio 
of the measurement can be increased by weighting 
each pixel by the optical depth. We find that the 
increase in signal to noise is approximately 10%. 
However, when filtering with the MF this weight- 
ing scheme increases the contamination from ther- 
mal SZ, since the thermal SZ signal is also propor- 
tional to the optical depth. We therefore do not 
weight by the optical depth to reduce the risk of 
contaminating our result. 

In principle the correlated errors on the three 
dipole to values can be calculated, with covariance 
matrix 

n = (x T wxy 1 x T wcw T x(x T wx)- 1 

(19) 

where C is the pixel-pixel noise covariance ma- 
trix. Since it is not practical to compute a ma- 



trix as large as C, we calculate the covariance us- 
ing simulations. The dominant source of error is 
CMB and instrument noise, which we estimate by 
passing 100 CMB and noise realizations through 
our pipeline, performing the least squares fit on 
each. There are additional sources of error from 
our uncertainty in the optical depth, and from the 
random component of the galaxy cluster peculiar 
velocities, which are calculated from equation 10. 
We estimate the errors from both of these terms 
by passing simulations of the kSZ signal through 
our analysis pipeline, performing the least squares 
fit on each realization. The scatter in the derived 
values of (3 then provides an estimate of the noise 
correlations between the dipole directions, which 
is then used to calculate the x 2 significance of the 
measured dipole in each redshift shell. We find 
that the errors from the latter two terms are neg- 
ligible compared to the CMB and noise error, and 
that the errors on the dipole directions are Gaus- 
sian with a narrower distribution in the direction 
perpendicular to the galactic plane. This is due 
to clusters lying within the galactic mask being 
removed from the fit. 

As a check on this process, we also estimate the 
errors from the WMAP maps by re-calculating the 
dipole fit after rotating our weights map away from 
the clusters. We rotate the weights map in incre- 
ments of 10' in the galactic longitudinal direction. 
We remove points within 1° of the cluster center 
to reduce any residual SZ emission. This leaves us 
with 2147 different directions. This process allows 
us to preserve the angular distribution of the clus- 
ter sample in the fit. When either filter is used 
we find errors that are consistent with those found 
from the CMB noise realizations. The errors we 
find by rotating the weights map are less reliable 
than the those calculated from the CMB and noise 
realizations, since when we rotate our weights map 
some regions get rotated into the galactic mask, 
which we then exclude from the fit. The dipole 
is therefore calculated from fewer clusters in this 
scheme, and so the errors we quote are from the 
CMB and noise realizations. 

We find that 100 CMB and noise realizations is 
sufficient to estimate our errors. Taking a typical 
error for one of the dipole coefficients in the Wiener 
filtered maps to be 0.8/iK we find that the standard 
deviation of the sample variance is 0.09/xK for 100 
noise realizations. Even if our errors are 10% too 
large we would still not detect a significant cluster 
dipole with our filters. 
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4.5. Conversion to Velocity Dipole 

We create a conversion matrix to calculate the 
velocity dipole from the temperature dipole, 



a v = Mai 



(20) 



where a v are the velocity monopole and dipole co- 
efficients, &t are the temperature coefficients and 
M is a 4 x 4 conversion matrix. We calculate M by 
creating four simulated maps of kSZ signal alone 
from our cluster sample. In one map all of the 
clusters are given a monopole velocity. Each of the 
other three maps has a bulk flow velocity in one 
of three basis directions. These maps are passed 
through our analysis pipeline to use the same fil- 
tering and velocity fitting processes that are used 
for the real data. The four fit coefficients from 
each of the input maps are the elements of each 
row of M _1 . The matrix M is close to diagonal 
in all redshift shells, and the velocity and tem- 
perature dipoles are close to alignment, as they 
should be. To convert from /iK to 1000 km/s we 
use the following values for the diagonal entries of 
M: 4.07 (our Wiener filter), 3.41 (KAKE Wiener 
filter), 0.075 (matched filter), and 0.12 (unbiased 
matched filter). For the dipole contribution: 



M = 



10 



-12 



TCMB T e ff 



1 
1 
1 



(21) 



where T c g is the effective optical depth after fil- 
tering and we have used the same normalization as 
KAKE (which makes M a factor of smaller 
than for the conventional definition of a\ m ). For 
our Wiener filter we find T off = 2.8 x 10" 5 . This 
value is smaller than the average optical depth of 
our cluster sample (r = 4.9 x 10~ 3 ) due to the fil- 
tering, as well as the fact that the signal is diluted 
by beam smoothing and averaging over the aper- 
ture used to calculate the cluster dipole. We find 
a lower effective optical depth than the value used 
in the KAKE analysis of T e g <~ 10" 4 , which could 
explain why we find velocity limits ~ 10 times 
larger. We are confident that we correctly recover 
the cluster velocities and show in figure 13 that we 
correctly recover the velocities in simulations. 

4.6. Tests of the Method 

We have performed tests to check that the re- 
sults from our filtering and dipole fitting pro- 
cedures are not contaminated by systematic ef- 
fects. First we verify that we can correctly re- 
cover a known cluster bulk flow from simulated 



kSZ maps alone, then we consider maps contain- 
ing kSZ, CMB and instrument noise. In section 5 
we describe the effects of thermal SZ, unresolved 
sources, and galactic emission on this process. 

Using simulated maps with a kSZ component 
only, but that include a cluster bulk flow we have 
verified that we recover the input bulk flow veloc- 
ity using our method. We assign all of the clusters 
in the simulated maps a kSZ signal from both the 
bulk flow and a random velocity drawn from the 
expected ACDM distribution in equation 10. For 
each choice of input bulk velocity we generate 25 
realizations, filter the maps, and fit for a cluster 
dipole using the same pipeline that is used for the 
WMAP data. We repeat the process with three 
different input bulk flow directions: (galactic lat- 
itude, longitude) = (0°, 0°), (0°, 90°) and (90°, 
0°) and bulk flow velocities of km/s, 500 km/s 
and 1000 km/s (7 cases with 25 realizations each). 
Figure 12 shows the results of the fits. We do 
not plot the recovered velocity for the cases where 
the simulated bulk velocity is non-zero because we 
use our kSZ simulations to calibrate the recovered 
velocity, and so the mean recovered velocity is ex- 
act by design. The scatter in the recovered dipole 
direction is due to the random component of the 
cluster line of sight velocity (with variance given 
by equation 10) and uncertainty in the cluster op- 
tical depth. The recovered amplitude in the maps 
with no bulk flow velocity has significantly larger 
scatter when it is calculated using the UF. This re- 
flects the lower signal to noise of this filter caused 
by the extra degree of freedom. 

We have repeated the test using simulations 
that additionally contain CMB and instrument 
noise. We create 100 realizations for each WMAP 
channel and add a bulk flow signal to the kSZ com- 
ponent, with a range of velocities logarithmically 
spaced between 100 and 100,000 km/s in the direc- 
tion of galactic latitude —11° and longitude 103°, 
which is the bulk flow direction found by KAKE 
in the redshift shell extending to z — 1. Figure 13 
shows the recovered bulk flow amplitude and fig- 
ure 14 shows the error in the direction of the re- 
covered dipole when all clusters are included in the 
fit. The lower signal to noise of the UF is apparent 
in both of these figures. Figure 15 shows the re- 
covered monopole velocity in a different set of 100 
simulated WMAP realizations. 



15 



Matched Filter 



Unbiased Filter 



' E 
< 




1000 



0.05 0.10 0.15 
Mean redshift 



0.05 0.10 0.15 0.20 
Mean redshift 



Fig. 12. — Left: Error in the recovered dipole di- 
rection from simulated kSZ maps with a bulk flow 
velocity of km/s (top), 500 km/s (middle) and 
1000 km/s (bottom) in the directions: (latitude, 
longitude) = (0°, 0°) black lines, (0°, 90°) blue 
lines and (90°, 0°) red lines, using the MF. Right: 
Same for the UF. The arrows are the 95% confi- 
dence limits. 
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Fig. 14. — Error in the direction of the bulk flow 
velocity in simulated maps containing CMB, noise 
and kSZ with the MF (black) and UF (blue), when 
all clusters are included in the fit. 
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Fig. 13. — Recovered bulk flow velocity in simu- 
lated maps containing CMB, noise and kSZ using 
the MF (black) and the UF (blue). All clusters 
are included in the fit. The x-axis is the bulk ve- 
locity input into the simulated maps, the y-axis is 
the recovered velocity. The red line indicates per- 
fect recovery. When the input bulk flow is small, 
a dipole is not detected and the data points pro- 
vide an estimate of the scatter in the recovered 
amplitude. 
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Fig. 15. — Error in the recovered monopole veloc- 
ity in simulated maps containing CMB, noise and 
kSZ with the MF (black) and UF (blue). 
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5. Systematic Effects 
5.1. Thermal SZ 

There are not expected to be intrinsic large scale 
moments in the tSZ signal, but due to the limited 
size of our cluster sample the tSZ dipole signal 
caused by random scatter is not negligible. We 
have examined the effects of thermal SZ on our 
results using our own simulations of the tSZ emis- 
sion from our cluster sample, which are described 
in section 4.2, as well as a simulated tSZ map pro- 
duced using the Planck Sky Model, described in 
section 3.2. 

We find that when we fit for the cluster dipole 
amplitude in our tSZ simulated maps filtered with 
the MF, using all clusters in our sample, there is 
a 3(7 bias in the y and z directions, which equates 
to a dipole amplitude of approximately 4000 ± 400 
km/s in the direction of l,b = (88°, 54°) ± 40°, 
which is 66° from the KAKE bulk flow direction 
in the same redshift shell. However, in the PSM 
tSZ map we find a reduced dipole of 2500 km/s, 
in the direction of I, b = (327°, 10°), which is 137° 
from the KAKE bulk flow direction. These simu- 
lations suggest that contamination to the cluster 
dipole measurement is at the la level and so can- 
not be ignored. As expected, we find that the tSZ 
monopole in the simulated maps is large and equiv- 
alent to a kSZ signal from clusters with velocities 
of 12000 km/s in our simulations, and 10000 km/s 
in the PSM maps. 

The simulated tSZ maps filtered with the UF 
have a cluster monopole and dipole amplitude that 
are more than an order of magnitude lower. In our 
simulated maps containing only tSZ we find a re- 
covered cluster dipole amplitude of approximately 
200 km/s. In the PSM simulated tSZ map we find 
a recovered cluster dipole amplitude of 70 km/s, 
significantly less than the measured cluster dipole 
and noise. The monopole amplitude is 200 km/s in 
our simulated tSZ map and 100 km/s in the PSM 
tSZ map. 

Although the UF removes the tSZ signal suf- 
ficiently for our purposes, it does not remove it 
to the ACDM limit. The residual tSZ in the fil- 
tered maps, and the cause of the non-zero clus- 
ter monopole and dipole is source confusion. The 
emission from clusters with small angular separa- 
tion overlaps more strongly in the Q band where 
the beam is large, than in the W band. When the 
channels are subtracted the tSZ signal does not 
subtract perfectly in the central cluster pixel and 
there is residual tSZ in the map. Since this only 



affects a small number of clusters, it leaves a dipole 
signal in the map. 

5.2. Radio Point Sources 

At frequencies below 100 GHz extragalactic ra- 
dio sources can be a significant contribution to 
CMB maps. We use simulations of the unre- 
solved extra-galactic radio point source emission 
produced by Colombo & Pierpaoli (2010) and sum- 
marized in section 3.2, to verify that our bulk 
flow measurement is not significantly affected by 
sources. Because the model is based on sources 
observed by NVSS, the simulated maps retain in- 
formation about the distribution of sources on the 
sky and account for the increased radio emission 
at galaxy cluster positions, caused by clustering of 
radio galaxies. 

The uncertainty on the radio point source 
monopole and dipole in the maps recovered using 
the UF are larger than in the maps filtered with 
the MF. In the maps filtered with the MF we find a 
non-zero dipole introduced by radio sources with 
a signal of amplitude 0.2a in each of the dipole 
directions. In the redshift shell encompassing all 
clusters we find signals of 0.2tr, O.lcx and — 0.3<7 in 
the x, y and z directions previously defined. We 
find a monopole signal of — 2a, equivalent to a kSZ 
monopole velocity of approximately —2000 km/s. 
In the maps filtered with the UF we find a larger 
signal with an amplitude of — 0.9c, la and 1. lo- 
in each of the dipole directions. The monopole 
signal is —6a, as expected it is large because the 
radio sources all add signal to the maps whereas 
the dipole is caused by sample variance. 

Because the UF is designed to remove tSZ by 
a specific weighting of channels, the radio source 
contribution is actually amplified by the filter. In 
the UF the Q band channel is given relatively more 
weight, so that emission with a tSZ spectrum is 
canceled when the channels are combined. The ra- 
dio point source signal is strongest in the Q band 
and weakest in the W band. The UF therefore has 
a large contribution from the Q band radio signal 
that is not offset by the W band signal, where the 
radio point source signal is much weaker. The re- 
sulting map generated using the UF therefore has a 
greater residual point source signal than the map 
generated with the MF. The radio point source 
signal and the tSZ signal could in principle both 
be suppressed simultaneously by using a filter de- 
signed to remove both signals (Tegmark & Efs- 
tathiou 1996, e.g.,). However, the increase in the 
number of degrees of freedom that this would in- 



17 



troduce would further reduce the signal to noise of 
the measurement. 



5.3. Galactic Emission 

For our analysis we use the WMAP foreground 
reduced maps which have been processed by the 
WMAP team to suppress the galactic signal out- 
side of the WMAP galactic mask (Jarosik et al. 
2010). We test the effect that any residual galac- 
tic foreground would have on our results by re- 
peating the analysis using maps that have not had 
the galactic components suppressed. The cluster 
monopole and dipole that we obtain from these 
maps is therefore an upper limit on any residual 
galactic contamination in our results. 

Figures 16 and 17 show the sum of the filtered 
PSM simulated maps. Both maps show residual 
galactic emission that is not visible in the WMAP 
maps (figures 7 and 8) due to their foreground re- 
duction procedure. The larger visible foreground 
signal in the PSM simulated map filtered with the 
UF is caused by galactic synchrotron. The syn- 
chrotron signal is smaller in the map filtered with 
the MF because the large scale galactic signal is 
removed at each frequency by the MF, but re- 
mains in the maps filtered by the UF. When the 
filtered maps at each frequency are combined, the 
signals from the galactic components are not can- 
celed since the filters are designed to only cancel a 
signal that has a thermal SZ spectrum. 

Because the maps containing only the individ- 
ual contribution of a particular foreground com- 
ponent are available for the PSM simulations, we 
have also run our analysis pipeline on each compo- 
nent separately to identify the source of any bias. 
The results when all clusters are included in the 
dipole fit are shown in Table 2. 

We find no significant bias in the monopole or 
dipole of the maps filtered by the MF. The com- 
bined map filtered by the UF has the largest galac- 
tic contribution from synchrotron emission, with a 
— 5.7ct signal in the monopole (labeled ' in the ta- 
ble) and a — 4tr signal in the dipole direction that 
points towards the center of the galaxy (labeled ^ 
in the table). The dust and free- free monopole sig- 
nals largely cancel resulting in a — 6.2a monopole 
signal. The significant dipole from synchrotron 
emission is partially canceled by the dust and free- 
free signals, resulting in a dipole signal that is not 
significant and would not be considered a detec- 
tion. 

As a further check, we repeat the cluster dipole 
analysis with the cluster coordinates rotated by 




Fig. 16. — Sum of the PSM simulated maps con- 
volved with the WMAP beams and filtered with 
the MF. This map shows galactic emission around 
the edges of the mask that is absent from the 
WMAP foreground reduced maps due to WMAP's 
foreground reduction method (Jarosik et al. 2010). 
The color scale is altered to better show the galac- 
tic emission by mimicking the effect of replacing 
the data by sinh _1 (data) (Gorski et al. 2005). 




Fig. 17. — Sum of the PSM simulated maps con- 
volved with the WMAP beams and filtered with 
the UF. The sinh -1 color scale that was used in 
figure 16 is used here. 
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one degree away from the true cluster locations. 
We find that the emission from the galactic compo- 
nents is not significantly changed. Since the emis- 
sion is not localized at galaxy cluster locations, 
we expect the galactic signals we find to be sup- 
pressed by the WMAP foreground reduction pro- 
cedure, and so we expect the galactic signals in 
Table 2 to be an upper limit on residual galactic 
emission in the WMAP maps. 

6. Results 
6.1. Dipole 

Table 3 and figure 18 show our KAKE filter 
pipeline results, table 4 and figure 19 show the 
WF results, table 5 and figure 20 show the MF re- 
sults and table 6 and figure 21 show the UF results. 
The points in figures 18-21 are from the WMAP 7 
year data. The green line is the noise bias, the 
red and blue lines are the 95% and 99.7% confi- 
dence limits that there is no bulk flow, which are 
estimated from our realizations containing CMB 
and instrumental noise. We find smaller errors for 
our Wiener filter than for the KAKE filter since 
our filter suppresses more power at I < 100 as 
shown in figure 24. The noise at these multipoles 
is dominated by the CMB and is correlated be- 
tween channels (Keisler 2009). We do not find a 
significant dipole in any of the redshift shells using 
any of the filters. In figure 21 some of the points 
are close to the 99.7% confidence limit. Since this 
filter is much less sensitive than the other filters 
which have results that are consistent with noise, 
this result cannot be due to a bulk flow. Instead 
it is likely that it is at least partly due to radio 
point source contamination which in section 5.2 
we estimated to cause a ~ la bias in each of the 
a\ m components in maps filtered by the UF. We 
have repeated our analysis for the WMAP 5 year 
data and find no significant cluster dipole. We find 
that for the MF the CMB and instrument noise 
contribute approximately equally. For the UF the 
noise is dominated by the instrument noise, with 
approximately 90% contribution. 

The UF is an order of magnitude less sensitive 
than the MF. This is due to the channels being 
combined in a way that is not optimal for maxi- 
mizing the kSZ signal, but results in cancellation 
of the tSZ signal. The reduction in sensitivity of 
the UF is largely due to the limited frequency cov- 
erage of the maps we use, and so an experiment 
with greater frequency coverage, such as Planck, 
would perform better with this filter. 

Although the dipole amplitude we find is consis- 
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Fig. 18. — Left: Cluster dipole amplitude in the 
maps filtered by the KAKE filter. Points with a 
plus sign have redshift shells with minimum red- 
shift of 0, those with a cross have shells with min- 
imum redshift of 0.05 and 0.12. The green line is 
the noise bias, the red line is the 95% confidence 
limit that there is no bulk flow and the blue line is 
the 99.7% confidence limit. Right: Cluster dipole 
amplitude in km/s. 
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Fig. 19. — Cluster dipole amplitude in the maps 
filtered by the WF. 
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Table 2 

monopole and dipole in the psm simulations.* 



V^UlllJJ) Wilt. 11 L 


Filter 


Monopolc 






Dipole b 






./7J7 T//K1 

V ^1 Lr txv J 


95% Confidence Limit \nK] 






oo [/iK] 


a a 


x [nK] 


<r a 


y [mK] 


a 


z [nK] 


a 






Synchrotron 


Matched 


0.40 


0.038 


0.054 


0.0020 


0.16 


0.0079 


-0.086 


-0.0048 


0.19 


62 


Dust 


Matched 


0.36 


0.034 


0.034 


0.0013 


0.38 


0.019 


0.018 


0.0010 


0.38 


62 


Free- free 


Matched 


0.042 


0.0040 


-0.46 


-0.017 


0.38 


0.019 


-0.042 


-0.0023 


0.60 


62 


All galactic 


Matched 


0.80 


0.077 


-0.37 


-0.014 


0.92 


0.046 


-0.11 


-0.0061 


1.00 


62 


Synchrotron 


tSZ Bias Removing 


-469 


-5.7+ 


-717 


-4.0" 


156 


1.0 


-120 


-1.0 


743 


432 


Dust 


tSZ Bias Removing 


58 


0.70 


170 


0.94 


-159 


-1.1 


70 


0.58 


244 


432 


Free- free 


tSZ Bias Removing 


-101 


-1.2 


188 


1.0 


177 


1.2 


-29 


-0.24 


260 


432 


All galactic 


tSZ Bias Removing 


-513 


-6.2 


-357 


-2.0 


174 


1.2 


-80 


-0.66 


405 


432 



*In the redshift 0.00-1.00 shell. 

a The uncertainty estimates are obtained from CMB and noise realizations. 

b Thc dipole basis directions arc x (I = 0°, 6 = 0°), y (I = 90°, b = 0°), z (I = 0°, b = 90°). 
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Fig. 20. — Cluster dipole amplitude in the maps 
filtered by the MF. 
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Fig. 21. — Cluster dipole amplitude in the maps 
filtered by the UF. 



tent with zero, our limits to the bulk flow velocity 
are tighter in some directions than in others. Fig- 
ure 22 shows the 95% confidence upper limit to 
the bulk flow over the whole sky using the results 
from the MF with all clusters included in the fit. 
Table 7 gives the upper limit to the flow in the 
direction of other well-known dipoles. The upper 
limit in all directions is above the limit expected 
by cosmic variance in the ACDM model, and is 
above the measured low redshift flow. 

6.2. Monopole and Higher Moments 

The monopole result from the MF is expected 
to be strongly contaminated with thermal SZ emis- 
sion, and we find a result consistent with tSZ con- 
tamination: —6868 ± 838 km/s. In the maps fil- 
tered with the UF we find a monopole consistent 
with zero: 4692 ±8947 km/s. We find no moments 
between multipoles 2 and 5 that are significantly 
different from zero. 

6.3. Comparison with SuZIE measure- 
ments 

The SuZIE II experiment has placed limits on 
the bulk flow using pointed observations of galaxy 
clusters (Benson et al. 2003; Benson 2004). SuZIE 
II made simultaneous measurements of the SZ ef- 
fect in three frequency bands, centered at 145 GHz, 
221 GHz, and 355 GHz, from the Caltech Submil- 
limeter Observatory, with a 1'.5 beam FWHM at 
all frequencies. A significant detection of the ther- 
mal SZ signal was made in 15 clusters; no signif- 
icant detection of the kinetic SZ signal was made 
in any of the clusters observed. 



Benson et al. (2003) fit a thermal and kinetic 
SZ spectrum to the measured cluster temperature 
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Fig. 22. — 95% confidence upper limit to the bulk 
flow in the redshift 0—1 shell. We find a bias in 
simulated thermal SZ maps equivalent to a bulk 
velocity of ~2500 km/s. 
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Table 7 

95% Confidence Limits on the Dipole. 



Dipole Name 


I [deg] 


b [deg] 


Velocity [km/s] 


Reference 1 


CMB 


263.99 


48.26 


3485 


[1] 


Lauer & Postman 


220 


-28 


4625 


[2] 


Watkins 


287 


8 


4509 


[3] 


KAKE 


267 


34 


3886 


[4] 



iReferences: [1] Jarosik et al. (2010), [2] Lauer & Postman (1994), 
[3] Watkins et al. (2009), [4] Kashlinsky ct al. (2009) 



was found in the SuZIE cluster sample. We find 
that the SuZIE bulk flow limits are more sensitive 
than our results with the MF. We have combined 
the cluster dipole from the SuZIE data with the 
dipole limit we find in the WMAP data. When 
combining the data we weight the dipole ampli- 
tude in each of the x, y, and z directions by the 
inverse noise variance. We find that the 95% con- 
fidence limit in the direction of the CMB dipole is 
increased from the SuZIE result of 1500 km/s to 
1800 km/s with the combined data, because our 
best fit dipole points in a different direction to the 
SuZIE dipole. With the UF our result is given 
less than 1% weight, and the combined result is 
unchanged from the SUZIE dipole measurement. 
The cluster monopole velocity in the SuZIE data 
using all of the clusters is —570 ± 320. The result 
with the UF is not sensitive enough to tighten this 
constraint. 

7. Conclusions 

We have used the kinetic Sunyaev-Zeldovich ef- 
fect to look for large scale moments in the galaxy 
cluster line of sight velocity distribution in the 
WMAP 7 year data. We use a multi-frequency 
matched filter that maximizes the cluster signal 
to noise ratio. We use a sample of 736 clusters 
derived from the ROSAT X-ray catalogs, and cal- 
culate the dipole at the locations of the clusters in 
the filtered maps. 

We find no evidence of a cluster dipole in the 
WMAP 7 year data in any of the redshift shells we 
use, consistent with predictions from the ACDM 
theory. Using kSZ simulations we create a matrix 
to convert the temperature dipole amplitude to a 
velocity amplitude. We find a 95% confidence up- 
per limit to the flow of 3485 km/s in the direction 
of the CMB dipole, and 3886 km/s in the direc- 
tion of the KAKE claimed flow (galactic longitude 
267°, latitude 34°). We have performed our anal- 



ysis on the WMAP 5 year data used in the KAKE 
analysis and find no evidence of a cluster bulk flow. 
We find that results obtained using our Wiener fil- 
ter have greater sensitivity than the results pro- 
duced using the KAKE filter, and our matched 
filter gives results approximately ~ 3 times more 
sensitive than those from the KAKE filter. The 
reason for the difference between our velocity lim- 
its and those of KAKE is due to the increased sen- 
sitivity of our pipeline and the different factor we 
use to convert from fiK to km/s, which depends 
on the kSZ simulations used. 

Our analysis differs from that of KAKE in sev- 
eral ways. We use a different filter to suppress the 
CMB component. We have tried three filters: a 
Wiener filter, a matched filter and a matched fil- 
ter that suppresses the tSZ emission. We use a 
different cluster sample. Our sample contains 736 
clusters outside of the galactic mask. The sam- 
ple described in KAKE has 674 clusters and the 
sample in Kashlinsky et al. (2010) contains 985 
clusters. For our Wiener filter we do not detect 
a cluster signal with any aperture we use. The 
numbers we quote in table 4 are from a 15' aper- 
ture which we found gives a higher signal to noise 
measurement than larger apertures. 

Using simulations of the tSZ signal we find that 
the results with our matched filter are contami- 
nated by thermal SZ, although the contamination 
is below the WMAP sensitivity. The ACDM model 
does not predict an intrinsic dipole in the tSZ emis- 
sion, but due to the relatively small size of our clus- 
ter sample, we find a non-zero signal. We estimate 
the signal to have an amplitude equivalent to a kSZ 
signal with a bulk flow velocity of — 2000 - 4000 
km/s. The KAKE analysis is performed with large 
~ 30' apertures around each cluster and so the tSZ 
signal is diluted, as seen by the small monopole 
values KAKE observe. 

When the maps are filtered with the MF we 
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find a monopole which is consistent with ther- 
mal SZ simulations, with a magnitude of (— 91.2 ± 
11.1) fjK, equivalent to a kSZ signal of -6868±838 
km/s. We use a modified multi-frequency matched 
filter that utilizes the different spectral shapes of 
the kSZ and tSZ signals to remove the thermal SZ 
bias. However, the signal to noise of the cluster 
dipole measurement is reduced by almost an order 
of magnitude, a consequence of using only three 
frequency bands in our analysis. 

The tSZ bias removing filter also has increased 
contamination from extra-galactic unresolved ra- 
dio emission. A filter could be constructed to sup- 
press this signal, however, since this would further 
reduce the signal to noise of the measurement, we 
do not further modify our filters to remove it. 

The limits we place on the cluster bulk veloc- 
ity can be decreased by using higher signal to noise 
measurements of the kinetic SZ signal, as well as by 
using a larger cluster catalog. We expect that our 
method can be applied to data from the upcoming 
Planck experiment. Planck is currently surveying 
the sky at a resolution of 5' at 217 GHz. In ad- 
dition Planck will itself produce a cluster catalog, 
which can be used for the analysis. Planck's wide 
frequency coverage from 30-857 GHz, including a 
channel at the thermal SZ null of 217 GHz, will 
allow the tSZ signal to be removed with a smaller 
impact on the kSZ signal to noise ratio. The ef- 
fect of radio point sources will likely be negligible 
because the strongest cluster signal is at frequen- 
cies where the radio point source signal is small. 
However, extra-galactic infrared sources with a ris- 
ing spectrum will probably be an important con- 
taminant. Our filter could be modified using spec- 
tral information about the IR source population to 
suppress emission with this frequency dependence. 
A paper describing the application of our method 
to the Planck experiment is in preparation (Mak 
et al. 2011). 
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Table 3 

Results from the KAKE filter.* 





Zmax 


{*) 


^median 


2<T 


N cl 


Monopole, ao [/iK] a 


aix [/iK] 


ai y \pK\ 


aiz [/iK] 


Dipole Amplitude, VUi [/iK] b 


0.0 


0.02 


0.014 


0.015 


0.0048 


28 


-1.4 ± 3.8 


-12. ± 10. 


-12. ± 8.7 


-6.6 ± 8.7 


18. 


0.0 


0.025 


0.016 


0.016 


0.0057 


36 


-6.2 ± 3.5 


-8.6 ± 8.7 


-11. ± 7.6 


-5.6 ± 5.5 


15. 


0.0 


0.03 


0.019 


0.019 


0.0075 


50 


-5.8 ± 3.1 


-1.8 ± 6.3 


-4.3 ± 6.3 


_2.4 ± 4.4 


5.2 


0.0 


0.04 


0.027 


0.030 


0.0100 


95 


-1.6 ± 2.2 


-3.2 ± 4.9 


-4.9 ± 4.3 


-5.6 ± 3.1 


8.1 


0.0 


0.05 


0.033 


0.035 


0.012 


139 


-2.3 ± 2.1 


-2.1 ± 4.9 


-6.0 ± 3.8 


-1.5 ± 3.0 


6.6 


0.0 


0.06 


0.039 


0.041 


0.015 


192 


-3.9 ± 1.9 


-2.4 ± 3.8 


-4.9 ± 3.5 


0.017 ± 2.8 


5.5 


0.0 


0.08 


0.050 


0.052 


0.019 


294 


-4.6 ± 1.7 


-1.2 ± 3.0 


-3.9 ± 3.0 


-1.9 ± 2.5 


4.5 


0.0 


0.12 


0.067 


0.066 


0.029 


445 


-4.4 ± 1.5 


-0.82 ± 2.5 


-2.7 ± 2.3 


-1.1 ± 2.2 


3.0 


0.0 


0.16 


0.080 


0.076 


0.039 


546 


-3.9 ± 1.5 


-0.77 ± 2.2 


-2.7 ± 2.1 


-0.86 ± 1.9 


3.0 


0.0 


0.20 


0.092 


0.083 


0.049 


619 


-3.6 ± 1.4 


0.024 ± 1.9 


-2.6 ± 2.0 


-0.31 ± 1.7 


2.6 


0.0 


1.0 


0.12 


0.097 


0.079 


736 


-3.6 ± 1.2 


0.55 ± 1.7 


-2.1 ± 1.8 


-0.26 ± 1.6 


2.2 


0.05 


0.30 


0.13 


0.12 


0.064 


578 


-3.9 ± 1.3 


1.2 ± 1.7 


-1.4 ± 1.9 


-0.098 ± 1.7 


1.9 


0.12 


0.30 


0.19 


0.18 


0.048 


271 


-2.1 ± 1.5 


2.9 ± 2.2 


-1.8 ± 2.4 


0.56 ± 2.0 


3.4 



Table 4 

Results from our Wiener filter.* 





Zmax 


(z) 


^median 


2 CT 


N cl 


Monopole, ao [/iK] a 


ai x [a*K] 


ai y [/iK] 


aiz [/iK] 


Dipole Amplitude, yfC{ [/iK] b 


0.0 


0.02 


0.014 


0.015 


0.0048 


28 


0.40 ± 2.5 


-4.8 ± 5.3 


0.73 ± 5.6 


-0.97 ± 5.0 


5.0 


0.0 


0.025 


0.016 


0.016 


0.0057 


36 


-1.9 ± 2.2 


-3.8 ± 4.7 


-2.1 ± 5.0 


-3.2 ± 3.7 


5.4 


0.0 


0.03 


0.019 


0.019 


0.0075 


50 


-3.0 ± 1.9 


0.066 ± 3.9 


0.15 ± 3.9 


-1.4 ± 2.8 


1.4 


0.0 


0.04 


0.027 


0.030 


0.0100 


95 


-2.5 ± 1.3 


-0.58 ± 2.3 


-1.7 ± 2.5 


-1.6 ± 1.9 


2.4 


0.0 


0.05 


0.033 


0.035 


0.012 


139 


-2.2 ± 1.0 


-0.92 ± 2.1 


-1.3 ± 2.0 


0.096 ± 1.5 


1.6 


0.0 


0.06 


0.039 


0.041 


0.015 


192 


-3.4 ± 0.91 


-2.3 ± 1.6 


-1.8 ± 1.7 


2.1 ± 1.3 


3.6 


0.0 


0.08 


0.050 


0.052 


0.019 


294 


-3.8 ± 0.79 


-2.3 ± 1.4 


-1.9 ± 1.5 


-0.59 ± 1.0 


3.0 


0.0 


0.12 


0.067 


0.066 


0.029 


445 


-3.5 ± 0.62 


-1.6 ± 1.1 


-1.1 ± 1.1 


-0.29 ± 0.78 


2.0 


0.0 


0.16 


0.080 


0.076 


0.039 


546 


-3.7 ± 0.55 


-1.3 ± 0.91 


-1.1 ± 0.96 


-0.41 ± 0.68 


1.7 


0.0 


0.20 


0.092 


0.083 


0.049 


619 


-3.7 ± 0.51 


-0.91 ± 0.85 


-0.66 ± 0.91 


-0.59 ± 0.63 


1.3 


0.0 


1.0 


0.12 


0.097 


0.079 


736 


-3.7 ± 0.45 


-0.87 ± 0.80 


-0.95 ± 0.82 


-0.47 ± 0.61 


1.4 


0.05 


0.30 


0.13 


0.12 


0.064 


578 


-4.1 ± 0.50 


-1.1 ± 0.96 


-0.93 ± 0.91 


-0.74 ± 0.68 


1.6 


0.12 


0.30 


0.19 


0.18 


0.048 


271 


-4.0 ± 0.66 


-0.19 ± 1.6 


-1.2 ± 1.4 


-0.99 ± 1.2 


1.5 



"The monopole and dipole arc calculated within 15' of the cluster centers. In tables 5 and 6 the monopole and dipole are calculated in the central pixel 
of the clusters. 



Table 5 

Results from the matched filter. 



Zmin Zmax (z) z median z a N ct Monopole, ao \pK] a ai x [fj,K] aiy [/iK] a\ z \jjK] Dipolc Amplitude, y/C\ \pK] b 



0.0 


0.02 


0.014 


0.015 


0.0048 


28 


59 ± 65 


-39 ± 117 


-10 ± 140 


-56 ± 101 


69 


0.0 


0.025 


0.016 


0.016 


0.0057 


36 


55 ± 57 


-111 ± 110 


-79 ± 126 


-68 ± 70 


152 


0.0 


0.03 


0.019 


0.019 


0.0075 


50 


-15 ± 44 


-82 ± 97 


-110 ± 101 


-1.8 ± 60 


137 


0.0 


0.04 


0.027 


0.030 


0.0100 


95 


-22 ± 33 


-20 ± 60 


-107 ± 70 


-16 ± 40 


110 


0.0 


0.05 


0.033 


0.035 


0.012 


139 


-41 ± 27 


-45 ± 58 


-44 ± 55 


13 ± 36 


65 


0.0 


0.06 


0.039 


0.041 


0.015 


192 


-62 ± 24 


-16 ± 49 


-23 ± 43 


33 ± 29 


43 


0.0 


0.08 


0.050 


0.052 


0.019 


294 


-63 ± 18 


-39 ± 37 


-38 ± 32 


-7.1 ± 24 


55 


0.0 


0.12 


0.067 


0.066 


0.029 


445 


-73 ± 13 


-22 ± 28 


-6.1 ± 28 


-0.059 ± 20 


23 


0.0 


0.16 


0.080 


0.076 


0.039 


546 


-81 ± 12 


-18 ± 28 


1.8 ± 25 


-1.4 ± 18 


19 


0.0 


0.20 


0.092 


0.083 


0.049 


619 


-85 ± 12 


-14 ± 26 


5.7 ± 23 


-9.3 ± 17 


18 


0.0 


1.0 


0.12 


0.097 


0.079 


736 


-91 ± 11 


-5.7 ± 24 


-1.3 ± 21 


-8.7 ± 17 


11 


0.05 


0.30 


0.13 


0.12 


0.064 


578 


-99 ± 12 


-2.2 ± 26 


9.5 ± 24 


-16 ± 19 


19 


0.12 


0.30 


0.19 


0.18 


0.048 


271 


-112 ± 17 


11 ± 37 


-0.13 ± 37 


-30 ± 27 


32 



Table 6 

Results from the tSZ bias removing filter. 





Zmax 


(z) 


^median 


Za 


N ci 


Monopole, ao [/uK] 


aix [m k ] 


ai y [/xK] 




Dipole Amplitude, y/C{ [;uK] b 


0.0 


0.02 


0.014 


0.015 


0.0048 


28 


118 ± 425 


-242 ± 869 


-1045 ± 866 


-29 ± 746 


1073 


0.0 


0.025 


0.016 


0.016 


0.0057 


36 


-156 ± 326 


269 ± 756 


-1268 ± 862 


-449 ± 537 


1373 


0.0 


0.03 


0.019 


0.019 


0.0075 


50 


73 ± 267 


1256 ± 598 


-29 ± 610 


-93 ± 420 


1259 


0.0 


0.04 


0.027 


0.030 


0.0100 


95 


-223 ± 193 


938 ± 410 


-281 ± 384 


159 ± 314 


992 


0.0 


0.05 


0.033 


0.035 


0.012 


139 


-164 ± 169 


616 ± 394 


-310 ± 298 


170 ± 241 


710 


0.0 


0.06 


0.039 


0.041 


0.015 


192 


-22 ± 140 


694 ± 344 


-200 ± 263 


-74 ± 224 


726 


0.0 


0.08 


0.050 


0.052 


0.019 


294 


-60 ± 122 


660 ± 280 


-230 ± 228 


15 ± 166 


699 


0.0 


0.12 


0.067 


0.066 


0.029 


445 


-38 ± 102 


607 ± 225 


-200 ± 180 


-131 ± 133 


653 


0.0 


0.16 


0.080 


0.076 


0.039 


546 


-56 ± 95 


591 ± 199 


-147 ± 159 


-182 ± 135 


636 


0.0 


0.20 


0.092 


0.083 


0.049 


619 


-14 ± 89 


549 ± 183 


-123 ± 154 


-214 ± 128 


603 


0.0 


1.0 


0.12 


0.097 


0.079 


736 


40 ± 77 


486 ± 165 


-68 ± 140 


-178 ± 119 


522 


0.05 


0.30 


0.13 


0.12 


0.064 


578 


98 ± 92 


367 ± 185 


66 ± 150 


-320 ± 134 


492 


0.12 


0.30 


0.19 


0.18 


0.048 


271 


189 ± 129 


230 ± 272 


343 ± 230 


-418 ± 184 


588 



a The monopole is contaminated by foregrounds with levels consistent with simulations. 

b The errors have a one sided distribution and so we use a 95% confidence limit and do not state error bars. 
c The modified KAKE filter is not normalized (see discussion in text). 



A. Expected Cluster Velocity Dipole 

We follow the derivation for the density distribution in Peebles (1973). The line of sight velocity field can 
be expanded as 



a lm = J drr 2 dn r (/)(r)Y e * m (r)v(r) • r 



(Al) 



where r is the comoving radial distance, <fr(r) is the comoving number density of objects in the sample, 
which we assume to be isotropic, v is the object peculiar velocity and Ye m (r) are the spherical harmonics. 
The peculiar velocity at a given wavenumber is related to the over-density by 



v(k)=*/fT J(k)- 



(A2) 



where / = (a/D) dD/da w O^ 55 , fi m is the matter density parameter, a is the scale factor, D is the 
growth function and Hq is the Hubble constant. The expansion in terms of spherical harmonics is given by 



^ _ t 1 T r1 \ d 3 k P(k) 

' (27t) 3 k 2 



\a em \ 2 } = f 2 Hl 



d z rcj)(r)Y tm {Y) k- f e 



-ik-r 



167T 2 f 2 H 2 j J d 3 r<b(r)Y ern (r)k-r ^{-if j t ,(kr)Y t , m ,(l)Y;, m ,(T) 

P'm' 



(A3) 



where P(k) is the matter power spectrum and ji(kr) are the spherical Bessel functions. Then 



Ci = (\a im \ 2 ) 



^fH 2 J dkP(k) (^J drr 2 ^(r)g e (k 



r) 



where 



(A4) 



gtikr) 



2£+l 



[tM(kr)-(l+l)j e+1 (kr)} 



(A5) 



The selection function, (p(r), is estimated by calculating the number density of clusters in our sample in 
radial bins, which we then interpolate to give a smooth function and normalize. 



B. Comparison with KAKE Results 

B.l. Filter Pipeline Comparison 

We have verified that we can reproduce the KAKE pipeline by comparing our filtered maps with the 
publicly available maps used for the Kashlinsky et al. (2011) analysis. We find that the mean difference 
between the W band maps passed through our pipeline and the filtered maps from the KAKE pipeline is 
— 1.5 x 10~ 10 toK indicating that our pipeline can reproduce their results. As a further check we calculate 
the spectrum of the filtered Wl channel map, and the absolute value of the difference between our filtered 
map spectra and the spectra from the KAKE maps in figure 23. In figure 24 we show the spectra of the 
WMAP maps passed through the two pipelines. The signal to noise is low at I < 100, which is why both 
filters suppress these multipoles. The KAKE filter is larger at I < 100 because the filter does not remove 
the fluctuations in the map that are due to cosmic variance. This leaves noise in the maps that is correlated 
between channels (Keislcr 2009). By further suppressing the map at £ < 100 our filter reduces the noise 
in our dipole measurement by a factor of <~ 2 while not significantly affecting the signal, as can be seen in 
figure 11. 

Since we model our clusters as point sources we can calculate the profile of the cluster signal in the filtered 
maps. In figure 25 we show the integrated signal out to a given radius in the Ql, VI, and Wl channels after 
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filtering with the Wiener filter. The signal is normalized so that if no filtering were applied the result would 
approach unity at large radius. The fact that the signal peaks at 10— 15' justifies our choice of a 15' aperture 
within which to calculate the monopole and dipole. 



10° 
10" 2 



W1 Map Spectrum 
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200 



400 
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Fig. 23. — Spectrum of the Wl channel map filtered by the KAKE filter (black) and the difference between 
the filtered W band map spectra from our KAKE filter pipeline and the spectra of the publicly available 
maps used in the Kashlinsky et al. (2011) analysis. 
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Fig. 24. — Filtered Wl channel map spectrum from the KAKE pipeline and our pipeline. 
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Fig. 25. — Integrated signal from a beam smoothed and filtered cluster in our simulated maps. 

B.2. Cluster Dipole Comparison 

In this section we remove clusters with bolometric luminosity Lx < 2 x 10 erg/s from our sample and 
use the WMAP KpO mask instead of the 7 year extended temperature analysis mask to better compare with 
the KAKE results. As before we remove the monopole and dipole of the map outside of the mask before 
fitting for the cluster monopole and dipole. We find that the sample of clusters satisfying the luminosity 
cut and the redshift cut of z < 0.16,0.2,0.25 is almost identical to the sample used in the KAKE analysis. 
There are some differences in the REFLEX, BCS and BCSe catalogs for clusters with luminosities close to 
Lx = 2 x 10 44 erg/s, and some differences for clusters in the CIZA catalog that are near to the WMAP 
galactic mask boundary. We show the differences we find in table 8. 



Zmax 


Number of clusters 
with L x > 2 x 10 44 erg/s 


Number of clusters KAKE 
include but we don't 


Number of clusters we 
include but KAKE do not 


0.16 


131 


13 


11 


0.20 


205 


14 


18 


0.25 


269 


14 


20 



Table 8: Number of clusters satisfying the luminosity and redshift cut in our sample and the KAKE sample. 



To check that these differences do not change the results by a large amount we have calculated the cluster 
dipole using the KAKE filter with both the KAKE cluster sample and our own sample. We find the results 
shown in table 9 when averaging over the Q, V and W bands. 
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Temperature \pK] 


Cluster sample 




Zmax 


m 


Cm 


X 




y 


<7y 


z 




Ours 


0.0 


0.16 


-0.74 


2.05 


1.44 


2.87 


-7.04 


2.98 


0.04 


2.68 


Ours 


0.0 


0.20 


0.02 


1.48 


3.04 


2.21 


-6.23 


2.52 


1.93 


1.97 


Ours 


0.0 


0.25 


0.04 


1.21 


1.56 


1.92 


-4.17 


2.14 


2.08 


1.72 


KAKE 


0.0 


0.16 


-1.07 


2.07 


1.51 


2.96 


-7.00 


2.86 


0.48 


2.72 


KAKE 


0.0 


0.20 


-0.22 


1.46 


3.06 


2.22 


-6.20 


2.45 


2.08 


2.02 


KAKE 


0.0 


0.25 


-0.20 


1.22 


1.68 


1.99 


-4.25 


2.12 


2.19 


1.76 



Table 9: Comparison of the cluster monopole and dipole in our cluster sample with the KAKE sample. 

All of the results have less than 3er significance with some of the a\ y components having ~ 2.5a significance. 
The results averaged over the W band channels only (but with the same filter and cluster samples used to 
produce table 9) are shown in table 10. 





Temperature \pK] 


Cluster sample 


Zmin 


Zmax 


m 


<Jm 


X 


<*x 


y 


Gy 


z 


Cz 


Ours 


0.0 


0.16 


-1.64 


2.11 


1.11 


3.12 


-8.63 


3.24 


0.01 


2.84 


Ours 


0.0 


0.20 


-0.41 


1.59 


2.71 


2.42 


-7.73 


2.67 


2.08 


2.04 


Ours 


0.0 


0.25 


-0.49 


1.30 


1.03 


2.06 


-5.53 


2.28 


2.10 


1.81 


KAKE 


0.0 


0.16 


-1.64 


2.13 


0.96 


3.20 


-8.15 


3.11 


0.54 


2.88 


KAKE 


0.0 


0.20 


-0.48 


1.55 


2.60 


2.43 


-7.45 


2.60 


2.28 


2.12 


KAKE 


0.0 


0.25 


-0.60 


1.30 


1.06 


2.13 


-5.43 


2.26 


2.25 


1.87 



Table 10: Same as table 9 for the W band channels only. 



These results have greater significance than the combined channel results. We still find nothing with 
more than 3<r significance, although the z = 0.0 — 0.20 shell a\ y value has 2.9cr significance. The result 
from Kashlinsky et al. (2011) is shown in table 11. 





Temperature [fjK] 




m 


X 


y 


z 


0.00 0.16 
0.00 0.16 
0.00 0.16 


-1.47 
-0.32 
-0.44 


1.20 
2.83 
1.30 


-8.26 
-7.58 
-5.57 


0.38 
2.13 
2.10 



Table 11: Result from (Kashlinsky et al. 2011). 



which is similar to our result above. These results are more significant than those obtained using the full 
cluster sample. This cannot be fully explained by the higher luminosity clusters having larger optical depths. 
In table 12 we show the average optical depth and number of clusters in the sample both with and without 
the luminosity cut. 





With L x cut 


No Lx cut 


Zmax 


Number of Clusters r avg x 10 3 


Number of Clusters 


T~avg X 10 a 


0.16 


131 


6.3 


573 


3.7 


0.20 


205 


6.6 


649 


4.1 


0.25 


269 


7.0 


713 


4.5 


All z 


327 


7.0 


771 


4.9 



Table 12: Mean optical depth of our cluster sample. 



Although the average optical depth is higher with the luminosity cut, the noise (which is roughly pro- 
portional to l/V-Wciusters) is a ^ so higher and so we expect a similar sensitivity both with and without the 
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luminosity cut. However, when we calculate the cluster dipole in maps filtered with our Wiener filter with 
our cluster sample we find much lower significance as shown in table 13. 





Temperature [ji K] 


Cluster sample z min z max 


m a m 


x a x 


y °y 


z a z 


Ours 0.00 0.16 
Ours 0.00 0.20 
Ours 0.00 0.25 


-1.58 0.45 
-1.01 0.40 
-1.15 0.34 


-0.21 0.90 
-0.08 0.68 
-0.84 0.57 


-1.75 1.09 
-1.05 0.76 
-0.98 0.68 


-0.35 0.58 
-0.51 0.48 
-0.62 0.49 



Table 13: Result from our filter pipeline with the luminosity cut. 



We conclude from these results that we can accurately reproduce the KAKE pipeline and do not detect 
a bulk flow with greater than 3tr significance in the WMAP maps. The significance of the results is lower 
when we use our Wiener filter and the full cluster sample than when using the KAKE filter and the cluster 
sample that has low luminosity clusters removed. By comparing with table 3 we see that the results with the 
smaller cluster sample are more significant than when those with the full cluster sample. However, since we 
would expect any bulk flow to be of similar significance with either of the cluster samples we conclude that 
there is no significant detection of a bulk flow. 
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